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Abstract 

Collapse models are modifications of quantum theory where the wave function is treated as 
physically real and the collapse of the wave function is a physical process. This appears to introduce 
a time reversal asymmetry into the dynamics of the wave function since the collapses affect only 
the future state. This paper challenges this conclusion, showing that in three different examples 
of time asymmetries associated with collapse models, if the physically real part of the model can 
be reduced to the locations in space and time about which collapses occur, then such a model 
works both forward and backward in time, in each case satisfying the Born rule. Despite the 
apparent asymmetry of the collapse process, these models in fact have time reversal symmetry. 
Any physically observed time asymmetries that arise in such models are due to the asymmetric 
imposition of initial or final time boundary conditions, rather than from an inherent asymmetry in 
the dynamical law. This is the standard explanation of time asymmetric behaviour resulting from 
time symmetric laws. 


* daniel. b edingham® philosophy, ox. ac. uk 
i owen.maroney@philosophy.ox.ac.uk 


1 




I. INTRODUCTION 


The idea of a physical wave function undergoing collapse seems to be at odds with time 
reversal symmetry. For example, during the collapse process an initially dispersed wave 
function might become spontaneously localised about some position. By contrast, the time 
reverse of this process would describe an initially localised state spontaneously dispersing. 
Collapse models treat the collapse of the wave function as a genuine physical process and 
therefore inherit this time reversal asymmetry. The most well known examples are the 
GRW model [1] and the CSL (continuous spontaneous localisation) model [21 [3] (for general 
reviews of collapse models see BlEl). 

There are grounds to be suspicious about the time asymmetry of wave function collapse 
since the Born rule can be understood to apply both forward and backward in time as 
outlined in the two-state vector formalism (TSVF) of Ref. [6]. The usual way in which 
quantum mechanics is applied involves the construction of ensembles of pre-selected states. 
This allows us to make predictions about the future state. The TSVF shows that by post 
selecting the hnal state we can also make retrodictions using the same laws. This implies 
that the basic laws of quantum mechanics can be understood in a time symmetric way and 
that time asymmetry results from the way in which we choose our statistical ensembles. 

In collapse models the physical collapse of the wave function is a manifestation of the 
Born rule. However, it is a manifestation of the Born rule only for the forward in time 
direction. Given that the Born rule can be stated in a time symmetric manner, we must be 
careful that we are not imposing our own prejudices about how the future differs from the 
past. 

Here we will argue that collapse models can indeed be understood in a time symmetric 
way. In order to do this we make use of a number of ideas: 

o Collapses happen randomly: A typical feature of collapse models is that the wave 
function undergoes spontaneous collapses in some preferred basis. This process may 
be discrete or continuous. The probability for any given collapse to happen should 
be given by the Born rule. For example, in the GRW model [1] the preferred basis 
is the particle position basis and the collapse process is discrete. The wave function 
'iptixi,... ,xn) for N distinguishable particles usually satisfies the Schrodinger equation 
but from time to time, and with hxed probability per unit time for each particle, it 
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makes a jump of the form 


A+ = ( 1 ) 

Each particle has its own random sequence of jumps occurring at a different set of 
random times. The jump operator j is of the form 

jU) = exp {-xy2a^) /(iTaY\ (2) 

with a some hxed length scale. The action of j is therefore to quasi project the wave 
function for particle i about some position 2 ; (the set {j‘^{z — Xi)\z G M} forms a 
POVM). The collapse centre ^ is chosen randomly from a probability distribution 

/ dxi - ■ ■ dxN\j{z - Xi)'ipt\^ 

/ dxi ■ ■ ■ dxN\'^t\^ 

This is precisely the Born rule probability for a quasi projection of the form j. The 
hxed rate of collapses can be chosen such that individual particles are rarely affected, 
but a bulk mass with large numbers of particles suffers frequent jumps. In this way 
macroscopic pointers rapidly commit to dehnite readings. 

o A model of wave function collapse can be given an ontology in terms of the locations of 
the collapses: This is usually referred to as the hash ontology [7]. The wave function 
is an object that does not exist in ordinary position space. However, the jumps in 
the GRW model are localised in space and time. It therefore makes sense to treat the 
collapses themselves, or the locations of the collapse centres z, as the basis for the 
local beables of the theory - the mathematical counterparts to real world events [8]. 
On a hne-grained scale the world appears as if composed from many discrete points. 
The local density of these points give a representation of the location of matter. In 
this picture the role of the wave function is in determining the probabilities for the 
various collapse locations. 

In some other collapse models, notably the CSL model and the lattice model discussed 
below the jumps act on a quantum held variable and occur throughout space 

and time. The resulting collapse centres take the form of a classical stochastic held. 
Again we expect to recover a picture of the world from this classical stochastic held 
by some coarse graining procedure. 
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o The wave function can be abandoned altogether fWf : As pointed out in Ref. HH, if 
a collapse model is considered as a stochastic law for generating the collapse centres, 
then the wave function can be relegated to the initial time from which it does not need 
to evolve. The collapsing evolution of the wave function corresponds to an updating, 
conditioned on the history of realised collapses, of the rule for determining the prob¬ 
ability of future collapses. In Ref.[TT] it is demonstrated that, given a set of collapse 
data generated by some specihc initial wave function, a generic wave function (subject 
to certain constraints on conserved quantities) will, after undergoing those same col¬ 
lapses, tend towards the same state as the original wave function. This implies that 
even the initial wave function can be replaced by a sufficient period of collapse data. 
The evolving wave function is then just a convenient calculation tool for making the 
theory Markovian. 


On the basis of these ideas, we conjecture that, given a valid set of collapse centres, 
we can form a picture of a collapsing wave function evolving both forward and backward 
in time, and in each case, the locations of subsequent collapses will satisfy the Born rule. 
In this sense a collapse model can have time reversal symmetry. The forward going and 
backward going pictures will not be the same for the same set of collapse data. In particular 
the backward going wave function will be affected by collapse events in the future and not 
the past. Nonetheless, the collapse centres, the local beables of the theory, will be consistent 
between the two pictures. 

To investigate this conjecture we will examine three different cases where collapse models 
might be thought to show distinctive time asymmetries: in Section the lattice model 


of collapse proposed by Dowker, Henson, and Herbauts P-ITT]: in Section hn the QMUPL 
(quantum mechanics with universal position localisation) model of Diosi [13] ; and in Section 
0 the generic tendency of collapse models to increase energy over time. In the hrst two 
cases we choose an initial condition, an initial wave function, and evolve forward in time in 
order to generate a set of collapses. We then show that there is an equivalent, backward 
in time picture of an evolving wave function undergoing the same set of collapses in the 
reverse order. We will show that the fixed collapses which were generated by the forward 
in time dynamics are distributed as though they had been generated by the backward in 
time collapse dynamics. In the final case we will show how the apparent monotonic increase 
in energy is compatible with time symmetric dynamical laws and discuss the signihcance of 
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asymmetric boundary conditions. We end with a summary in Section |Vj 


II. LATTICE MODEL 

Here we outline the lattice model of collapse proposed by Dowker & Henson [9] and 
further investigated by Dowker & Herbauts [101 m, in a slightly different presentation. 
We demonstrate some pertinent features of the lattice model. In particular, given only 
the stochastic field of collapse data, it is not possible to observe that the state is in a 
superposition of different preferred basis states. To discern the presence of some matter 
in a region, we must coarse grain to remove background noise. It turns out the scales on 
which we need to coarse grain are also the scales on which a superposition state collapses. 
This is important for our time symmetric picture since we should not be able to distinguish 
superposition states which appear in one time direction but not the other. We present the 
results of a statistical test, designed to demonstrate whether the collapses are distributed as 
though generated by the backward in time collapse dynamics. 

The model is a modihcation of the massive Thirring model on a (1 + 1)D null lattice 
[13]. This is a unitary fermionic quantum field theory in 2D Minkowski space. The lattice 
is shown in Figjlj 



FIG. 1: The (1 + 1)D null lattice showing a space-like surface a, a vertex v, and a link /. 

In the time direction (vertical) the lattice extends arbitrarily far into the past or future. 
In the spatial dimension there are N vertices denoted by black circles and we impose periodic 
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boundary conditions so that the space-time is in fact a cylinder. The centres of the links are 
also indicated by black squares. These are the locations of the collapse events. A quantum 
state I To-) is dehned with reference to a space-like surface a. This surface must cut through 
the links of the lattice as shown in the hgure. The surface must also satisfy the periodic 
boundary conditions. 

We can divide the lattice into columns representing the spatial location which we label 
i = 1,2,, 2N as shown in Figjl] As we pass along a space-like surface a we cut through 
each of the 2N columns once. To each column we assign a qubit with basis states |0), 
corresponding to empty, and |1), corresponding to occupied; the full quantum state is a 
tensor product of each of the qubit states so that the basis states of the system are of the 
form \ui,U 2 ,... ,U 2 n) where the Ui G {0,1} represent the individual qubit basis states. For 
example, a state of the form |1, 0,1,1, 0, 0, 0,1) on a lattice of size A^ = 4 can be understood 
either as a particular quantum held conhguration or a state of 4 ‘particles’ at 4 specihc 
locations. The vacuum, or unexcited quantum held state is |0, 0, 0,0, 0, 0, 0, 0). We will 
refer to this basis as the preferred basis since it is the basis into which the system tends to 
collapse. The Hilbert space for the quantum state has dimension 2^^. 

The model also makes use of a classical stochastic held a taking values «/ = 0 or 1 at 
random on each link I of the lattice. 



FIG. 2: Elementary evolution I: crossing a vertex. 

There are two types of elementary evolution of the state which occur when the surface a 
advances in one of two diherent elementary ways. One corresponds the surface a crossing a 
vertex. This is shown in Figj^as the surface passes over the vertex v from a to ex'. If the 
vertex involved in this elementary evolution is one which connects links at positions i and 
i + 1 (modulo 2N), we dehne a unitary operator by 

U.^ = ll (g) • • • 0 li_i g) Ui^i+i ® li+2 g) • • • g) l2Af, (4) 
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where Ui^i+i is some 4D unitary operator acting on the ith and {i + l)th qubits. In crossing 
the vertex v the state changes according to 

(5) 


n 



FIG. 3: Elementary evolution II; crossing the centre of a link. 

The second type of elementary evolution occurs when the surface passes over the centre 
of a link. This is shown in Figj^ where the surface passes the centre of the link I from a' to 
a". If the link is located at position i on the lattice we dehne a jump operator by 

Si{ai) = li ® • • • 0 li-i 0 Ji{ai) 0 Ij+i 0 • • • 0 l2N- (6) 

In crossing the link I the state changes according to 

1^.") = (7) 

The jump operator describes the collapse of the wave function. Without the jump operator 
this model is the light-cone lattice massive Thirring model of Ref.[Tl]. Depending on the 
value of ai, the operator Jj acting on the ith qubit takes the form 

M0)= G^ {|0.)(0.|+A1l,)(li|} 

^.(1) = {-’f I0.X0.I + |1.){1.|} (8) 

with X the fixed collapse parameter, 0 < X <1. The jump operators satisfy 

(0) + JKi) = 1. (9) 

These are Krauss operators. They form a POVM on the ith qubit state. If X is close to 
0 then whenever they act they effectively perform a projective measurement on the qubit 
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state in the 0,1 basis. If X is close to 1 their action nndges the qnbit state slightly towards 
either the 0 or 1 state with collapse reqniring many snch jump operations. 

The value of the held variable controls whether the collapse favours the |0) state or 
the |1) state. This held variable is chosen randomly and the probability that the held takes 
value ai on the link I is given by 


cr = 


( 10 ) 


This rule ensures that the jump operators act in accordance with their Born rule proba¬ 


bilities. This is because the probability (10) is precisely the Born probability for the quasi 
projection Jj. 

From these elementary rules we can derive the rules for evolution from a general state on 
an initial surface Uj to a hnal surface <7/ in the future. We suppose that in getting from cTj 
to (7/ we must cross n links and m vertices. The links we label {/i,..., ln\ and the vertices 
we label {ui,..., Um}- It follows that the hnal state is of the form 




( 11 ) 


where T is the time ordering operator. The probability for the held values ..., ai^} is 


given from (10) by 


1 ) ■ • • 5 ^In 1 ^*) 


( 12 ) 


{'I',. I®,,) 

This probability is well dehned given only the partial ordering of links and vertices imposed 
by the space-time causal order [9]. 


h a unitary operator of the form 

m 
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i sin 6 
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0 

0 


0 \ 
0 
0 

1 / 


(13) 


where /F'\ denotes an incoming state where both the ith and the {i + l)th qubits are in 
the |1) state; ^ ■ denotes an incoming state where the fth qubit is in the |1) state and the 
{i + l)th qubit is in the |0) state; etc. Here, 9 controls the speed of the ‘particles’ on the 









lattice with, e.g., 9 = 0 corresponding to light speed particles and 9 = 7i/2 corresponding 
to stationary particles. All excitations of the quantum held travel with the same speed 
although they may move either to the left or the right. This unitary operator conserves the 
number of particles (the number of occupied qubits). 

If we eliminate the jumps from the dynamics we can unwind the state evolution backward 
in time using the inverse unitary matrices 

l't.,.) = r[uy...u„y]l>i-„,), (i4) 

where T is the anti time ordering operator. Taking the complex conjugate gives 

= (15) 

demonstrating that the complex conjugate state obeys the usual unitary dynamics backward 
in time. 

Now we include the jumps and suppose that the backward in time dynamics satisfy 


|4;.) = r ■ ■ ■ U.„| !$• ). 


(16) 


This looks the same as the forward in time rule 0 other than that the time ordering has 
been replaced by anti time ordering. In general |*ho-i) will not be the same as Itko-i). Our 
aim is to show that the field values {cx/j,..., ai^} (which are fixed by the forward evolution) 
also satisfy the probability rule 


1 ,..., aijaf) - 


/<j)* \ 

\ (Tj I CTi / 

/\T/* |\T/* \’ 

\ af\ af/ 


(17) 


the counterpart to Eq.(12). If this is the case then the backward in time evolution uses 
precisely the same dynamical rule as the forward in time evolution. The wave functions for 
the two cases will differ in general on any given surface a but the stochastic held a, the basis 
for local beables, will be consistent. 


A. Coarse graining and collapse time scale 

Given a set of collapses (the held a in our lattice model), the forward in time wave 
function undergoing these collapses is likely to look (at least at the micro level) diherent 
from the backward in time wave function. If the time symmetric picture of physical collapse 
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is to work, it must be the case that we are not able to observe the state of the wave function 
from the collapse data in a way which would allow us to distinguish these two cases. 

In this section we adapt a calculation given in Ref. HD] and show that the time scale 
necessary to observe the presence of an excited quantum held on the lattice is the same time 
scale on which a superposition of such excited states will collapse. This means that it is not 
possible to directly observe a superposition in the preferred basis given only a. 

Suppose that the system is in the vacuum state |0, 0, ■ ■ ■ ,0). The evolution of the system 
dehned by Eqs (|^ and ([^ will not change the state. The stochastic held a will take value 
0 with probability 1/(1 + and 1 with probability X^/(l + X^). We assume that X has 
a value close to 1. This is natural since we do not want single particle states to rapidly 
collapse. Consider a region of space-time R containing M links. Within this region the 
mean and variance of the held value are given by the binomial distribution 

1 + X2 ’ ^ M(l+X2)2- 

Now suppose that there is some non-vacuum state with mean held value aji in R. If we are 
to be able to observe against the background noise we require 

O' <C \aR — /i|. (19) 

If we write e = 1 — X, taking e to be small, we hnd that \aR — fi\ is at most (9(e) (occurring 
when the state is maximally excited in the region R). Since a ~ we therefore must 

have 


M > e“l (20) 

This determines the size of the region necessary to be able to identify an excited state. 

Now consider a block of n qubits in the 1 state with all other qubits in the 0 state, 
superposed with a disjoint block of n qubits in the 1 state with all other qubits in the 0 
state. For simplicity we assume that 9 = 7r/2 so that the unitary evolution has no ehect. 
Write the state as |R) -|- \B). 

After we have evolved for a number m of time units the state will be of the unnormalised 
form X^^\A) + X^^\B), where Ma + Mr = 2nm = M, and Ma/b is the number of links 
for which the held ai corresponds to the A/B state qubit eigenvalue on that link wherever 
the A and B state qubit eigenvalues are different. 
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Since the value of e is small, the field takes values 0 or 1 with probability ~ 1/2. The 
division of M into Ma and Mb has a distribution with standard deviation in Ma/b oc y/M. 
This means that one of the states will be suppressed with respect to the other by a factor 
of ~ exp(—e\/M). Therefore there is exponential suppression when 

M ~ e“l (21) 


We now put this together with result (20). If the ontology of the lattice collapse model is 
given by the stochastic held a then any feature observable against background noise requires 
coarse graining over a region containing ^ e~^ links. Since a superposition state is only able 


to survive for 


,-2 


links then we cannot directly observe the huctuating densities during 


the collapse process. 

Consider the example above and suppose that the state collapses to the \B) state after 
a certain period of time. If we then evolve the state backward in time it will simply stay 
in the \B) state. This can be consistent since the period during which the state was in a 
superposition of |y4) and \B) when viewed forward in time is brief enough that it cannot be 
distinguished from \B) if we only have access to the stochastic held a. 


B. Time reversed collapse on the lattice 


To test for consistency between forward in time and backward in time collapse dynamics 
on the lattice we apply the following test: 

1. We generate a some held data by starting with some initial state |To) on an initial 
time slice t = 0 and evolving forward in time. This results in a random held a. 


2. We then reverse the evolution in time using Eq.(16) applying the collapses again in 


the reverse order using the same held data a. For each link I we make note of the 
probability that the held ai takes the value 1 conditional on all a to the future of 
/. These probabilities can take values anywhere in the range [0,1]. Once we have 
returned to time f = 0 we have a set of reverse time probabilities for the held on each 
link on the lattice to take value 1 alongside the realised set of held values a. 

3. We divide the set of reverse time probabilities into a set of bins with boundaries 
0 < Pi < p 2 < • • • < 1- For each bin j we take the average probability as pj = 
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{pj-i — Pj)/ 2 and we count the number of collapse events in that bin, rrij. Using the 
binomial distribution we determine the mean number of times that we expect to see 
the stochastic held take value 1, pj, and the variance in this number, (t|. We use 
a simple test to check that the data should be approximately normally distributed 
(essentially that mj should be sufficiently large given the average probability pj for 
the bin) and discard those bins where this is not the case. We then count the actual 
number of times Uj that the held a for each of those events realised the value 1. From 
these we calculate a chi-squared statistic 


X 


2 


E 


i>h 



( 22 ) 


4. Using the theoretical distribution for the chi-squared statistic we calculate a p-value 
giving the probability that the held data or something more extreme could have been 
generated by the reverse time probabilities. In statistics a p-value of less than 0.5% is 
a typical standard for ruling out a hypothetical model. 

Figure]^ shows a typical example of a simulation with 2N = 16 in which we set X = 0.5 
and 6 = 7t/4. The initial condition is of the form 


^'o) 


10 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 0 ). 


(23) 


This is a single ‘particle’ state for the quantum held. We evolve for 100 time steps. We 
choose an arbitrary total ordering of links and vertices such that elementary evolutions sweep 
from left to right. (To be precise we start on a given time slice, we then evolve past the 
leftmost vertex, followed by the link above it to the left and then the link above it to the 
right. We then do the same for the vertex to the right and so on until we have evolved 
the time slice by one time step. The procedure is then repeated.) Each pixel in the hgure 
represents a link on the lattice so that each graph is 16 pixels wide and 100 pixels high. 

The left hand panel shows the expectation value (To-|A;|To-) where 

A; = li (g) • • • (g) li_i g Aj g Ij+I g • • • g l 2 N- (24) 


and 


A = |U)(U 


(25) 
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FIG. 4: (a) Quantum field expectation value for the forward in time state, (b) Stochastic 
field a. (c) Quantum field expectation value for the backward in time state. 


This is the expectation of the eigenvalue of the zth qubit. The surface a is the surface 
immediately following evolution across the link I given the total ordering which we use. We 
use a grey scale to plot the values of (\l/o-|A/|To-) with black corresponding to the value 1 and 
white corresponding to 0. A clear particle trajectory winding around the lattice is apparent. 

The central panel shows the realised held values a on the lattice generated randomly using 
the probabilities calculated from the forward evolving state IT^). Here black corresponds 
to a; = 1 and white corresponds to ai = 0. Here the particle trajectory is less apparent. 
As discussed in Sec jll A| a coarse graining procedure is necessary to eliminate background 
noise and establish the particle’s whereabouts to within space time regions of greater than 
1 pixel. 

The right hand panel again shows (\l/o-|Ai|To-), but this time the state evolves backward 
in time using the reversed total ordering and with collapses generated by the same held a. 
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The particle trajectory is very similar to that in panel (a). 

A careful observation of FigQa) reveals that dispersion occurs in the upward direction: 
localised particle states tend to become diffuse with increase in t before recollapsing with 
dead ends fading away. The same happens in the opposite time direction in Figj^c). Despite 
these differences the overall picture of a particle moving about on the lattice is broadly 
consistent. 



FIG. 5: Distribution of p-values for 500 simulations 


As outlined above we calculate the chi-squared statistic (22) and a p-value for the held 


data in panel (b) under the null hypothesis that data was generated by the backward in 
time collapse dynamics. In Figj^we plot a histogram of the p-values generated by running 
the simulation 500 times. The distribution of p-values is approximately uniform over the 
range [0,1]. This is as we would expect if the null hypothesis is correct. 

We conclude that held data generated by forward in time collapse dynamics is distributed 
as if it had been generated by an equivalent backward in time collapse dynamics. In this 
sense the collapse dynamics satishes a form of time reversal symmetry. This suggests that 
a backward in time moving observer could use the same mathematical procedure, based on 
a collapsing wave function to determine the past state of the stochastic held, as a forward 
in time moving observer would use to estimate the future state of the held. 
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III. TIME REVERSAL SYMMETRY FOR A LOCALISED WAVE PACKET 


We now turn our attention to the QMUPL model [13] for a particle in a localised state. 
When the particle is localised the collapse dynamics reduce to a classical diffusion process. 
We demonstrate that this process is asymmetric with respect to time reversal, but show it is 
still possible to construct an equivalent, backward in time process resulting from the same 
collapses occurring in the reverse order. 

The QMUPL model can be thought of as a continuous version of GRW describing wave 
function collapse for distinguishable particles. It is also a limit of the CSL (continuous 
spontaneous localisation) model when the particle density is low and the collapse length 
scale is large compared to the length scale of a wave packet [15] . In the QMUPL model the 
state vector for a single particle satisfies a quantum state diffusion of the form 

) = ^-iHdt - y(x - {x)tfdt + |(x - (x)t)d5t| iV't). 




Here Bt is a standard Brownian motion and the collapse parameter g controls the rate at 
which collapse of the wave function occurs. 

A typical feature of continuous collapse models is that after a sufficient period of time, 
the wave packet of an individual isolated particle achieves a stable localised shape. This 
happens when the dispersive effects of quantum theory balance with the localising effects of 
the collapses. When the particle is in this condition the dynamics are considerably simplified. 
The state can be characterised by only the central values of position {x)t = = Xt 

and momentum (p)* = {'ipt\'p\'^t) = Pt for the wave packet. The collapse dynamics causes 
these phase space parameters to undergo a classical diffusion process. For the QMUPL 
model this process is given by [T61 - IT8] 


dxt = —dt H— —dBt, 
m ^/m 

d-Pt = 


(27) 

(28) 


These equations describe a classical diffusion through phase space. In what follows we 
restrict our attention to this simplified model of a diffusing wave packet. This avoids much 
of the complication of an arbitrary wave function undergoing collapse. 


We first show that Eqs (27) and (28) do not have time reversal symmetry. The state of 
the system at any point in time is described by x, p, and B. Consider a sequence of two 
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states Si and S '2 at times ti = t and t 2 = t + At respectively. We write 


(29) 

(30) 


Si = {x,p,B} 

S '2 = {x + Ax, p + Ap, B + AB} 


Consider the change as we go from state Si to S 2 . The change in x is Ax, the change in p 


is Ap and the change in B is AB. From Eqs (27) and (28) we expect that 

1 


Ax = —At + 
m 

Ap = |a5. 


m 


AB, 


(31) 

(32) 


We next dehne a time reversal transformation T. For a classical phase space trajectory this 
involves a change in the sign of p simply becanse playing the movie of these events backward 
in time, the particle appears to move in the opposite direction to when the movie is played 
forward in time. The time reversed states are therefore given by 


Sl = {x^ iP^ iB'^} = {x,-p,B} (33) 

Sj = { 2 : + Ax, —p — Ap, B + AB} (34) 


Consider now the change in going from S 2 to Sf. The change in x, Ax^ = —Ax, the 
change in p, Ap"’" = Ap, and the change in B, AB^ = —AB. Inserting into Eqs (31) and 
(32) resnlts in 


Ax^ = —At + 


m 


AB"’ 


m 


Ap^ = -^AB^. 
2 


(35) 

(36) 


This is a different pair of eqnations from those which describe evolntion forward in time. 
This would appear to mean that it would be possible to determine the forward direction 
of time from a given sequence of states. For example, we could simply observe Ax and 
Ap and see if there is a positive or negative correlation in the spontaneous jumps. If the 
correlation is positive then the evolution is forward in time; if the correlation is negative 
then the evolution is backward in time. 

As with the lattice model we aim to show that this is just an artefact of the mathematical 
formulation and that it is possible to understand what is going on in a time symmetric way. 
Note that the expectation values x* and pt are features of the wave function. If we regard 
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the wave function as no more than a convenient way to encode the collapse history then it 
may be that we are adding in the time asymmetry by this procedure. Let us return to the 
idea that the collapse locations are fundamental. 

The reason for the stochastic motion of Xt and pt is that collapses are occurring randomly 
on either side of the centre of the wave packet causing it to spontaneously jump about. By 
taking the continuous limit of the GRW model we can show that the locations of the collapse 
centres are given by 


lAB, 

Zt = Xt - 7 —. 

9 At 

The collapse centres have a white noise distribution about the expected position xt- 

Now suppose that we only have information about Zt, i.e. the collapse locations form the 
basis of our local beables. We would like to conhrm whether it is possible to give different 
but consistent pictures of the evolutions of Xt and pt in either time direction, each satisfying 
the same dynamical law. We perform the following test analogous to the test carried out in 
SecJim 

1. We generate some collapse centre data Zi, i = 0,...,n by starting with an initial 
localised state characterised by the central position xq and momentum pq of the wave 
packet and evolving forward in time for n discrete time steps using the rule 



Xi+i = Xi + —At H— ]=ABi, 
m \/m 


Pi+i = Pi + ^^Bi, 

lABi 

Zi = Xi -;—. 

9 At 

Brownian increments ABi are generated randomly. This results in a 
phase space trajectory Xi, Pi. 


(38) 

(39) 

(40) 

forward in time 


2. Next consider the collapses from stage 1, Zi, in reverse order. Denote the backward in 
time phase space trajectory by x', p'. Starting at xjj = Xn and = —pn we use the 
dynamical law 


AB'_i = gAt{zi-i-x[), 

x'i_i = x[ + ^At + 

m ^/m 

p'_i = P- + |A5 '_i. 


(41) 

(42) 

(43) 
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In the first equation we back out the Brownian increments from the collapse centres 
using (37). The result is a backward in time phase space trajectory x', p'. 


3. We perform a Kolmogorov-Smirnov test on the set of implied increments ^B[/ A/At to 
see if they £t a normal distribution. This results in p-value for the reverse time AB' 
data. This tests whether the collapses occur with white noise distribution about x'. If 
they do then the backward in time trajectory satisfies the same dynamical law based 
on the collapsing wave function as the forward in time trajectory. 



FIG. 6: Phase space trajectories forward in time (black line), forward in time with p —)■ —p 

(dotted line), and backward in time (grey line). 


Figure shows an example of a trajectory though phase space using g = 20, m = 1 and 


At = 0.001. The solid black line is the forward in time trajectory generated using Eqs (31) 


and (32). The trajectory starts at time t = 0 at the black circle at a: = 0, p = 0 and ends at 


time t = 1 at the black square. There is a clear positive correlation in the stochastic jumps 
in X and p. The dotted line is the straightforward time reverse of this trajectory obtained by 
transforming p ——p. It is characteristically different from the forward in time trajectory. 
The correlation between stochastic jumps in x and p is negative. The grey line shows the 
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reverse time trajectory determined by the procedure outlined above. This trajectory starts 
at the grey square at time t = 1 and ends at the grey circle at time t = 0. This trajectory 
approximates the dotted line but with positively correlated jumps in x and p. Visually, there 
is no way to distinguish the micro dynamics of the backward trajectory from the forward 
trajectory. 



FIG. 7: Expected position with time. Black line is evolved forward in time; grey line is 

evolved backward in time. 

Figure is from the same simulated trajectory. It shows how x changes with time. The 
dots are the collapse centres generated by the forward in time dynamics and used again in 
the backward in time evolution. The black line is the forward in time expectation Xt- It is 
seen that the collapse centres are distributed about this line and both the collapse centres 
and Xt follow the same trend. The grey line is the backward in time x[ which is slightly 
different from the forward in time Xt but which also sits well within the distribution of 
collapses. 

Figure 1^ shows a histogram of p-values generated by running the simulation 5000 times. 
We hnd that the distribution of p-values is reasonably uniform in the range [0,1] indicating 
that indeed the increments IS.B' j belong to a normal distribution. 
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FIG. 8: Distribution of p-values for 5000 simulations. 

We conclude that the collapse data (the dots in Fig is distributed as though it had 
been generated by the backward in time collapse dynamics. 

IV. EQUILIBRIUM AND BOUNDARY CONDITIONS 


The two examples considered so far consider specihc models of objective collapse, exhibit¬ 
ing properties that look like an inherent time asymmetry (collapse events affecting the wave 
function after, but not before, the event; and positive correlation between the diffusion in 
mean position and momentum). Our third example looks at a more generic feature of col¬ 
lapse models: non-conservation of mean energy. This appears to give a directly observable 
time asymmetry. In doing so, we will shed light upon the role of initial and hnal boundary 
conditions in causing time asymmetric behaviour. 

It is well known (see, e.g. HI) that one of the physical effects predicted by collapse models 
is a gradual increase in energy. A monotonic mean energy increase can be seen to follow 


simply from Eq.(28). The diffusion in momentum of an individual particle is such that 


momentum is just as likely to go up by some amount as it is to go down by the same 
amount. However, since energy is a convex function of momentum, E ~ the energy 
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increases on average. It is important to understand how this can be the case with a time 
symmetric picture since this increase in mean energy only appears in one time direction. 

Consider a set of free localised particles, initially each with p = 0 at time t = 0. The total 
energy is zero. After a period of diffusion, the particles will have a Gaussian distribution of 
momenta centred on zero, but with a non-zero total energy. Viewed as a process backward in 
time, the energetic particles would appear to simultaneously tend towards zero momentum 
and energy. Even though we have shown that the collapse events are compatible with the 
Born rule applied to a backwards evolving wave function, the backwards in time evolution 
appears unlikely. For example, we would hnd that the collapses undergone by particles 
with positive momentum would tend to cause them to reduce momentum, and the collapses 
undergone by particles with negative momentum would tend to cause them to increase 
momentum. By conditioning on the sign of the momentum, a reverse time observer would 
appear to see a conspiracy, that leads to the mean energy decreasing. By observing the 
change in mean energy over time, the direction of time could be deduced. 

The conspiratorial behaviour for the reverse time picture comes about as a result, not 
of an inherent asymmetry in the dynamics, but as a result of the time asymmetric use 
of boundary conditions. From the forward in time picture, we have imposed the initial 
condition at time t = 0 that p = 0, and then tracked the behaviour for times t > 0. The 
evolution of the system then proceeds according to the Born rule, with no future boundary 
condition imposed. Eventually this will result in a diffuse nearly uniform distribution. 

Viewed in the reverse direction, we start with a diffuse distribution, evolving apparently 
according to the Born rule. We might expect this distribution to stay diffuse (or possibly 
spread out more if it is not in equilibrium). However, the behaviour we are tracking has 
a special feature that makes this impossible: these are pre-selected trajectories that must 
have come from localised p = 0 states at time t = 0. This condition on the trajectories 
introduces a pre-selective bias in the statistics of the reverse time collapse events, and leads 
to the apparent conspiratorial behaviour. 

To understand this more clearly, we will now address a well known problem for any 
stochastic theory [19]: such a theory cannot exhibit time independent behaviour (such as that 
required by the Born rule) in both temporal directions, unless the probability distribution 
is in equilibrium. 

The argument goes as follows: suppose a system follows a stochastic rule i?tptg(S'j|S'i) 
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for how it changes state from one time to another. At time to, the probability for state Si 
is PtQ^Si). At some later time ti > to, the stochastic evolution leads to the probability for 
state Sj of 

= (44) 

i 

From Bayes’ theorem, the system being in state Sj at time ti can be used to make retrod- 
ictions about the possible state at time to'. 

P„l„{S,\S,) = (45) 

If the future directed stochastic evolution lA*) is time independent (as we would 

expect from the Born rule), then adding any constant shift in time of r will not affect the 
transition probabilities: 


P,.+.|,.+.(Sj|S.) = ft.H,(S,|S.). 

However, imposing the same condition on the retrodictive probabilities 


(46) 


P.«+r|<.+x(S.|S,) = P,„|,.(Si|Sy) 


(47) 


gives 

P,.|,.(S,|Sj) = fl.,|.,(5,|5.) p‘°^t*^‘! 

^h+T[Oj) 

which can only hold if 

P,(5.) ^ P.(5;) ^ . 

Pt,.+AS,) P,,+r(S,) 


(48) 


(49) 


Normalisation of probabilities yields /(r) = 1. So the time independence condition cannot 
also hold in the backward direction, unless the system is in equilibrium with Pt^lSi) = 
PtQ+r{Si) = PE^Si). This would give a backwards in time retrodictive rule: 


Po\,,iS,\S,) 


Kit.iSAs,) 


Pe{S,) 

PE(s,y 


(50) 


which is time independant. 

We argue that this does not imply any sort of time asymmetry in the dynamics. Instead 
it is a result of the time asymmetric use of pre-selection statistics. Let us instead start with 
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the equilibrium distribution PE^Si), and then impose a pre-selection at time t = 0, that 
the system is in Sq. For times tf > 0, we have, as before, the forward in time predictions 
Rtf\o{Sj\So). If we now try to retrodict from the occurrence of Sj at time tf, to some 
intermediate time 0 < ti < tf, we get: 


Pt,\tEo{S^\SJ,So) = 


Rt,ltASj\Si)RtMS^\So) 


(51) 


In general, this will not lead to the reverse time retrodictive rule Eq.(50), nor will it be time 
independent. Pre-selecting the state Sq at time f = 0, biases the statistics of the retrodictive 
inferences in the time period 0 < < fj, as they must lead towards the state Sq as ti —)■ 0. 

However, exactly the same is true if we look at the forward in time statistics in the period 
tp < f_i < 0. First, if we impose the condition at time t = 0, that the system is in Sq, then 
at times fp < 0 we have the retrodictive inference: 


a.|o(Si|So) = (52) 

This leads to a time independent rule for tp < 0- By contrast, if we try to make predictive 
inferences in the period tp < t_i < 0, based on the occurrence of state Sj at time tp, we get 


Pt.i\tpfl{Si\Sj, So) 


Roit.AOiSi)Rt.,itp(SilS,) 

E,flo|i-.(0|Si.)fl,_.|,,(Sj.|S,) 


(53) 


Again, in general, this will not lead to the rule Rt_iitp(SijSj) nor will it be time independent. 
Post-selecting the state at time f = 0, biases the forward in time statistics in the period 
tp < f_i < 0, as they must lead towards the state Sq as t_i —)■ 0. 

There should be nothing terribly mysterious about the fact that an apparently time 
independent, future directed stochastic evolution law still does not yield time 

independent statistics when a post-selection is imposed at some time t > p. The post¬ 
selection biases the statistics. Equally so, therefore, there is nothing terribly mysterious 
that, when we start with a pre-selected ensemble, and look at the reverse time direction, the 
pre-selection introduces a bias in the retrodictive inferences. The apparent time asymmetry 
associated with the bias in the statistics is a result of the pre-selection, rather than anything 
fundamental to the dynamics. 

We may now return to the problem of non-conservation of energy in collapse models. By 
pre-selecting a set of free localised particles, with p = 0 at time f = 0, as an initial condition, 
we are biasing the sample for the reverse time statistics when f > 0. An analogous situation 
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can easily be constructed by starting with a wide and uniform distribution of initial momenta 
and post-selecting, at some time tf > t, only those trajectories which end up with momentum 
equal to zero at t = tf. The biases will now be present in the forward in time direction. 
Post-selected particles with positive momentum would tend experience collapses which cause 
them to reduce momentum, and those with negative momentum would tend to increase. 

Ideally, if the system has an equilibrium distribution, then to remove any biases we should 
start with samples whose distributions are approximately in equilibrium (see Ref. [20] h 
In the case of the QMUPL model, the momentum is not bounded, and an equilibrium 
distribution of momentum will have an unbounded mean energy, even if the mean momentum 
is hnite. Any selection at t = 0, with a hnite mean energy, will be followed in the forward 
time direction by a diffusion towards equilibrium that increases the mean energy for f > 0. 
However, when looking at times f < 0, the selection at t = 0 is a post-selection. The mean 
energy will be seen to be decreasing in time, converging on the value hxed by the post¬ 
selection. Once again, we see that the apparent asymmetry results from the use of initial or 
final boundary conditions, and not from any inherent asymmetry in the dynamics. A time 
symmetric interpretation of the stochastic dynamical law is therefore not necessarily at odds 
with an observation of energy increase due to collapse. 


V. SUMMARY 

We have shown that collapse models can be viewed in a way in which they exhibit a 
time reversal symmetry. This is perhaps surprising given the apparent time-directedness of 
a collapse of the wave function. The key is to treat the collapse centres or flashes as the 
fundamental stuff of the theory and the wave function as part of the dynamical law used to 
determine where the next flash will occur. Indeed we have argued that for a given set of 
flashes there are two equivalent pictures of a collapsing wave function: one going forward in 
time and one going backward in time, each collapsing according to the Born rule. 

This idea brings a new perspective to the problems of quantum theory. For example, in 
the usual picture of measurement the preceding quantum state is replaced by an eigenvector 
of the observable which is measured. Consider a particle passing through two misaligned 
polarisers. After passing through the hrst polariser the particle collapses to the hrst polari¬ 
sation state. It stays in this state for the duration of its trajectory before it passes through 
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the second polariser. In the backward in time picture the particle has the polarisation state 
of the second polariser at those points in time when it is situated between the two measure¬ 
ments. The resolution is to understand the wave function as an object which is part of the 
dynamical laws rather than part of the ontology. If we insist on understanding events in 
terms of the wave function then we must recognise that the wave function could be just as 
well determined by its future interactions as by its past m 
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